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1 Introduction 



One of the most debated issues in pre-big bang (PBB) cosmology concerns its initial condi- 
tions ^. How fine-tuned should these be in order to generate a Universe resembling ours? 
The principle of asymptotic past triviality (APT) |^ amounts to saying that, far enough into 
the past, the Universe was well described by a generic, perturbative solution of the low-energy, 
tree-level string effective action. In other words, the very early Universe is assumed to lie deeply 
inside the low-curvature, small-coupling regime. Under this assumption, for a critical super- 
string theory admitting the trivial Minkowski vacuum order by order in perturbation theory, 
the early Universe can be represented as a random superposition (a chaotic sea) of massless 
waves, propagating in all directions and with all (sub-string scale) frequencies. The particle 
content of the waves should represent all the massless degrees of freedom of superstring the- 
ory, i.e. each one of the possible marginal deformations of the two-dimensional CFT in trivial 
space-time. This can be hardly called a fine-tuned initial state! 

Understanding the evolution of such a Universe is clearly not a simple task. However, 
we know that, as long as the tree- level low-energy approximation remains valid, the classical 
field dynamics is both scale- and dilaton-shift-invariant since the string and Planck scales 
simply sit as overall factors in front of the action. General arguments suggest that, while most 
of the time these waves will just propagate linearly and independently, occasionally, through 
positive interference, overdense regions will form and, thanks to gravitational instability, will 
lead to gravitational collapse. It was argued that the interiors of sufficiently large and 
weakly coupled collapsing regions could give birth to Universes that could resemble our own. 
The above-mentioned classical symmetries guarantee that the scale of collapse will itself be a 
stochastic variable, whose distribution will be related to that of the original distribution of wave 
lengths. Since such a distribution naturally contains arbitrarily long (though not arbitrarily 
short) wavelengths, it looks very likely that large enough black holes would form with non- 
vanishing probability. 

Notwithstanding the appeal of these general arguments, it would be very instructive to see 
them at work in some explicitly soluble model. In this was partly done, in D = 4 for the 
spherically symmetric case, through use of some powerful results obtained by Christodoulou @] 
over more than a decade. More recently, exact analytic solutions were constructed by Feinstein, 
Kunze and Vazquez-Mozo (hereafter FKV) 0, who replaced the chaotic sea of waves by two 
colliding, homogeneous, planar-fronted waves. In the waves have infinite fronts and thus 
always lead to collapse on space-time scales L that are inversely related to the energy density 
(i.e. energy per unit of transverse area) in the waves: 

L ~ (G^piG^P2)-'/' . (1) 

While FKV only deal with D = A, we will show that the above result holds for any D. In the 
following we shall denote hjd = D — 2 the number of transverse coordinates. 

Although the case of infinite fronts is obviously an idealization, causality arguments lead 
to the conclusion that collapse takes place even for finite-front waves, provided the typical 



1 



transverse extension R of the waves is larger than L, i.e. if: 

i? > L ~ {GnPiG^p^)-'^' = i?" {GnE.GnE^)-'^' = = (2) 

i.e. provided Rs = {Gn\/s)~ > R, as naively expected. 

Other limitations of the FKV paper are that it deals only with gravitational and dilatonic 
waves and that, as we mentioned, it is restricted to D = 4. It has been pointed out recently 
that the approach to the cosmological singularity (which appears as the r = singularity 
inside the collapsing region) depends strongly both on the dimensionality of space and on the 
presence of other fields, in particular of the various p-forms that string theory possesses. It 
was argued, in particular, that turning on all the forms present in any consistent string (or 
M) theory, changes the monotonic Kasner-like behaviour of the gravi-dilaton system into an 
ever-oscillating behaviour a la BKL |^. 

As a first step in the direction of overcoming the two limitations of the work by FKV, we 
shall extend it both to an arbitrary number of (non-compact) dimensions and to the presence 
of the Kalb-Ramond B^i, field. Keeping the exact planar symmetry allows (see section 2) for 
an explicitly 0((i, (i)-invariant formulation of the dynamics in the string frame. While several 
properties of the solutions can be discussed in the general case (section 3), we have only been 
able, so far, to construct explicit analytic solutions in the case of 1?^,^ = and arbitrary 
dimensions (sections 4, 5), or when the 5^,^ field is generated through 0((i, (i)-transformations 
(section 6). This will suffice, however, to address some of the issues raised in Ref. [^. 



2 The 0{d, (i)-invariant field equations 

As explained in the introduction, we start with the low-energy, tree-level effective action of 
critical superstring theory in the string frame. Up to a classically irrelevant overall factor 

S = J dx^'^e-^ {r + d,ct>d^ct> - ^H^^pH^^'P^ , (3) 

where (p is the dilaton, R is the curvature scalar for the metric tensor ^f^j,, and Hap^ is the field 
strength for the antisymmetric field 

H^up = d^Bi,p + dpBpy + duBpp. (4) 

The problem we wish to study is the head-on collision of two infinite plane-front waves, 
which, without loss of generality, will be taken to move along the x^ axis. Assuming symmetry 
under translations along the d = D~2 spatial directions orthogonal to x^ (x* with i = 2, . . . D — 
1), it is obvious that the problem is endowed with d abelian isometrics. In a convenient 
coordinate frame, there will be no dependence upon the "transverse" coordinates x*, and thus, 
according to general arguments |10[, we expect to have an exact 0{d,d) invariance of the 
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classical field equations. The symmetry of the field equations (not to be confused with a true 
quantum symmetry) is most easily exhibited in terms of the invariance of a "reduced action" 
living in the non-trivial (here two-dimensional) subspace. In the case of 0{D — 1,D — 1)- 
symmetry (homogeneous Bianchi-I cosmologies) this was done in |1^ and led to a continuous 



extension of scale-factor duality of Bianchi-I cosmologies [|T2 



For the present purposes, we will adapt to our case a general result by Maharana and 
Schwarz |]T3, and write the reduced action coming from (Rf) as: 



S = dx^dx^y/^e~ 



R + do^Sd'^d) + ^Tr f doMT^d^M 



(5) 



whose notation we shall now explain. 



The equations of motion allow the metric g^^p to be taken block-diagonal, with blocks given 
by Qij and gap where roman indices («, j, . . .) will span the components of the tensors from 2 
to d+1 while the indices a and (3 take the values and 1. The explicit metric and curvatures 
appearing in (^ only refer to the latter. 

We arrange the components g^j in a dimensional matrix G. The matrix B will contain 
the components Bij of the antisymmetric field while the remaining components of B are set to 
0. M is then the 2 dimensional matrix defined by: 



M 



BG-^ 



-G-^B 
G-BG-^B 



Finally, the shifted dilaton is defined by: 



4> = 4> — - log det G + constant , 

where the constant will be conveniently fixed later. 

The reduced action (|^) is manifestly invariant under the transformations 

M -> n^MVL, 
where Vt denotes a global 0{d,d) transformation 

with 7] the 0{d, d) metric in off-diagonal form 



(6) 



(7) 



V 



Id 
Id 



M itself belongs to 0{d, d) and, being symmetric, satisfies 

M-qM = 7] , i.e. M"^ = r]Mr]. 



(9) 
(10) 



(11) 



(12) 



(13) 
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The above equations make the check of 0{d, (i)-invariance triviaL 

Let us now derive the (manifestly 0((i, (i)-invariant) field equations from (^. Varying the 
action with respect to the shifted dilaton gives 

R + 2(?"^Z}«D^0 - + ^^?°^Tr (dJA-^dfM) = 0, (14) 

while varying with respect to the 2-metric gap provides 

Rap + DaDp4> + ^Tr (dJA-^dpU) = 0. (15) 

Combining the trace of the latter equation with the previous one gives a simple equation 
for (f) which will play an important role later: 

L)"L)„exp(-0) = 0. (16) 

The variation with respect to M must be performed carefully because of the constraints 
coming from its definition. Following ||^, we take them into account by writing 

5M = n^MVL - M, (17) 

with VL = 1 + e m 0{d, d) . Expanding to first order in e, we obtain 

da [e-^V^g^^M-^dpM) = 0. (18) 

Equations. (|T^), ([T5|), (|T6D , (|18D compose the full set we wish to solve. Before doing so, let 
us write them more explicitly by going to the conformal gauge for the 2-metric 

9a(3 = e^r]ai3 (19) 
and by working with "light-cone" coordinates: 

Equations (|1^ , ([T5|), (p^, (|18D then simply become: 

9„9^exp(-0) = 0, i.e. dudycj) = duCpdyCp ; (21) 
(9„ [e-^M-^d^M) + dy [e-'^M-^d^M) = ; (22) 

dl4^ - duFdJp + -Tr (duM-^duM) = 0, same with u ^ v; (23) 
8 ^ ^ 

dud,'^ - dud.F + ^Tr [duM-'d,M) = 0. (24) 



Note that the two equations (psf ) have the form of Virasoro constraints of a two-dimensional 
CFT, while the other equations are of the evolution type. Furthermore, a straightforward, 
though not completely trivial, calculation shows that the integrability condition for eqs. (p3D 
holds and that eq. (P^ ) is just a consequence of the previous ones. Our discussion will thus be 
based on solving the set of equations (pl|)-(p3|). 
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3 General properties of the solutions 



In this section we will derive some general features of the solutions, which will be useful for the 
more detailed investigations to be described later on. 

The two colliding waves are defined to have their fronts at m = and v = 0, respectively, 
and thus to collide at u = v = (i.e. at x° = = 0). The two waves are not assumed 
to be impulsive, i.e. their energy density can have any (finite?) support at positive u and v, 
respectively. Space-time is thus naturally divided in four regions: 

Region I, defined by m, f < 0, is the space-time in front of the waves before any interaction 
takes place. It is trivial Minkowski space-time: 

ds'j = -2dudv + J2 {dx'Y ,B = , = 0o, (25) 

with a constant perturbative dilaton (exp(0o) ^ !)• It will be convenient to fix the constant 
in (^ as —00 so that, in region I, = 0. 

Region II, defined by m > 0, f < 0, is the wave coming from the left before the interaction. 
Metric and dilaton depend only on u and therefore the field equations allow us to take F = 
and write the ansatz: 



ds'jj 



2dudv + {u)dx^dx^ , Bi 



(26) 



Similarly, region III, defined by m < 0, f > 0, represents the wave coming from the right before 
the interaction. There, 



ds]ii = -2dudv + G"\v)dx'dx^ , B, 



kill/ 



(27) 



Note that we have not assumed any special shape for G, so that the results we give in this 
section will hold whatever the (relative) polarization of the waves. 



Finally, region IV (m > 0, f > 0) is the interaction region, with 



d^^ 



-2e^dudv + GlY dx^dx-' , 



(2^ 



and F, G^^ , (p^^ and B^^ are all functions of both u and v. Of course, the metric must be 
continuous along with its derivative on the boundary lines between the four regions. The same 
must be true for the dilaton and the antisymmetric field B. 

Let us begin by solving the equations in region II (and thus, by trivial analogy, in region III). 
The only non-trivial equation is the "Virasoro constraint", eq. (PS]), which (after momentarily 
dropping the subscripts II) reads 



It. 



(m^^m 



(29) 



where the dot indicates the derivative with respect to u. The r.h.s. can be written as the sum 
of three terms: 



1 



Tr 



(m-^m 



1 

4d 



Tr (G-^G) 



2 1 

+ 4 



Tr 



G-^G 



Tr 



G-^B 



(30) 
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where (^G ^G^^ is the traceless part of (^G ^G^. The first trace on the r.h.s. can be expressed 



in terms of the dilaton and the shifted dilaton as 

d 



1 

Ad 



Tr G-^G 




log det G 



(31) 



It is now useful to change variable from u to it, with 



Then, eq. 



becomes 



1 



'2 



d_ 

du 



1 

Ad 



-24>/d 



d_ 

du 



Tr 



G-^G 



Tr 



G^^B' 



(32) 



(33) 



where the prime denotes the derivative with respect to u. 



It can be easly proved that all terms on the r.h.s. of eq. ( p3D are negative definite, if 
we remember that G is symmetric with positive eigenvalues and B is antisymmetric. Hence, 
for any non-trivial wave, e"*^/*^, which is constant and identically equal to 1 in region_I, must 
acquire a non-vanishing, negative, and never increasing derivative in region II. Thus, e-^'"^ must 
vanish at some finite ib = ib* . Returning now to the coordinate m, we see that, if the dilaton is 
bounded (a necessary assumption if we want to use the tree- level effective action), there exists 
a finite u = u* where e"'^/'^ vanishes. Correspondingly, also det G vanishes, and the metric 
of the transverse space will collapse to zero proper volume, thereby producing a (coordinate) 
singularity. 

It is not too difficult to estimate the order of magnitude of u* by multiplying eq. ( p3| ) by 
g-</>/d g^j^^ integrating it once after having gone back to the original -u-variable. The result 
is an estimate of u*~^ and is given in terms of an integral over u of the energy density per 
unit volume, i.e. in terms of the energy density per unit area. One thus recovers the estimate 
u* ~ {Gnp)~^, in which the appearance of Newton's constant is somewhat fictitious since, for 
waves of a given geometry, the energy density scales like Gj^^. The final result, reported in 
eq. (|1|), is a Lorentz-boost-invariant way of writing the same expression once both waves are 
considered simultaneously. 



The same arguments can be repeated in region III, where det G has to vanish at some finite 
V = V* with (f) — > +00. These results generalize to any D and to non-trivial antisymmetric 
fields, a well-known result in D = A. 

Let us finally analyse region IV, where the interaction between the two waves occurs. We 
drop the subscript IV from all functions. We begin by using eq. (^Tj), which tells us that e~'^ is 
the sum of a function of u and a function of v. The unique function of this type that matches 
the boundary conditions with region I is 



(34) 
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We see that e must vanish on a hypersurface joining the coordinate singularities in regions 
II and III and contained in the region u < Uq, v < Vq within region IV. 

Let us now introduce two new sets of coordinates that simplify the analysis in region IV. 
One set is of the light-cone type: 

r = r{v) = 2e-^"^('') - 1, s = s{u) = 2e-^"(") - 1, (35) 

while the second set is of the t — x kind: 

e = ^(r + s)=e-^("''')~-t, (36) 
z = ^{s-r)= e-^"(") - e-^"^(^). (37) 

Note that the coordinates r, s run from +1 to —1 in region IV with their sum always positive 
except on the singular boundary where r + s = 0. Going from the original coordinates to either 
of the new sets changes only the conformal factor of the 2-metric. We may thus write, for 
instance, 

dsjy = -ef + dz^ + djdx'dx^ , (38) 
where / and G are functions of C, and z. In r, s coordinates, this becomes 

dsjy = —2e^drds + Gijdx^dx^ , (39) 

where / and G are now functions of r and s. The shifted dilaton is simply 

0= -log^ = -log(r + s)/2. (40) 



Finally, in these coordinates, eq. (^) becomes trivial and the only equations to be solved. 



1,1^31), become 



dr ((r + s)M-^dsM) + ds ((r + s)M-^drM) = 0, (41) 
(r + s)-^ + (r + sy^drF + ^Tr (dM'^dM) = 0, (42) 



(r + s)-^ + (r + s)-^dsF + -Tt (dM'^dM) = 0. (43) 

8 ^ ^ 

To end this section, let us discuss the relation between the above equations and those dis- 
cussed in Refs. p, where general solutions in the homogeneous and quasi- homogeneous 
[|T^] case were derived. Here we are dealing with fields depending also on one spatial coordinate. 
For this reason, we have built the matrix M by considering only the remaining d spatial di- 
mensions, excluding the gn component. In the approach to the singularity, we expect the time 
derivatives to dominate over the spatial ones (this hypothesis can also be tested and verified a 



posteriori, see section 6). Thus, close to the singularity, the general solutions of p|, should 
be recovered. To see the explicit connection between the equations of motion (^^3|) and those 
in P, O], we have to complete the M matrix with the missing rows and columns, extending 
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it to a 2 ((i + l)-dimensional matrix. We must also add the 5^11 component in the determinant 
of the metric used to construct the shifted dilaton and go over to cosmic time dT = e^/^d^. 
After some simple algebra, the equations of motion can be brought to a form identical to that 
of Refs. 1^, |I^, strongly suggesting that the asymptotic solutions of [|T^ will emerge near the 
singularity. 



4 Solutions with B = 



The case of i? = with parallel polarized waves was discussed in Ref. [^] for d = 2. We shall 
obtain below a generalization of their results to arbitrary d. For B = 0, M reduces to 



M 



G-^ 
G 



(44) 



In each region, we choose a diagonal G with Ga = e'^~^^\ where I^V'i = and A = ^logdetG 
goes to — oo at M = M* in region II and at v = v* in region III. 



The equations of motion follow easily from eqs. (|4lD and (^3]). They read 

1 



A,rs + 



2{r + s] 
1 



2 (r + s) 
1 r + s 



(A,, + A,,) =0, 

{'4'i,r + i^i,s) = 0, 



V + s) 4 
1 r + s 



r + s) 



^ {dxl + J2^ljl) = o. 



(45) 

(46) 
(47) 
(48) 



We see that the equations for the ipi and that for A are decoupled and can be solved 
separately. They are also formally the same, so the solutions, found by Szekeres 0, have the 
same structure: 



-(r + s)2ipi (r, s) 



ds 



p_l 

,s' 2 



1 + 2 



+ / dr' Ul + r')2tlj,{r',l) 



2 



1 + 2 



(1 -r) (s' - s) 
(1 + s') (r + s) 

1 — s) fr' — r 



(1 +r') (r + s) 



(49) 



The same expression holds for A, with the obvious replacements. In the above expressions, 
P_i{x) are Legendre functions written in standard notation. Once A is given, can be obtained 
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from eq. (|^, since is known. Finally, the function / (r, s) is given by an integral along a 
curve joining the point (r = s = 1), where it vanishes, to the generic point (r, s): 

{r,s) 

fir,s)= J [Irdr + lds,] (50) 
(1,1) 

with f r and f^s given by the r.h.s. of eqs. 
previously determined (p and ipi. 

We see that A and the ipi are singular on the hypersurface ^ = r + s = 0. So, in a very 
general way, we find that the collision of two plane waves leads to a (curvature) singularity in 
the space-time whatever the number of dimensions. Although we do not have, at the moment, 
such an explicit solution for the general case, the discussion given in the previous section makes 
us believe that a curvature singularity will always emerge along the hypersurface ^ = 0. 



47) and (Eq) as functions of (r, s) through the 



5 Asymptotic approach to the singularity 



For the solutions (^9]), the asymptotic behaviour for ^ — is easily found by taking the large- 
argument limit of the Legendre function [|] (see also Yurtsever 



i^i ~ e, (z) log^, 

~ k(^) log^, 
f{^,z) r^a{z) log^. 



(51) 
(52) 
(53) 



The coefficients multiplying the logarithm are functions of z, whose range, on the singular 
surface, is —1 — hi. One easily finds 



1 



ttVI + z 



ds 



+ 



1 



vrVl — z 



dr 



(l + s)^V^(l,s) 
(1 + r)^ (r, 1) 



s + 1\ 2 



s — z 



+ 



r + 1\ 2 
Vr + z, 



a {z) 



(54) 



(55) 



with K (z) given by the same expression as with ipi replaced by A. The sum of the must be 
zero according to the definition of the ipi. 



The asymptotic form of the metric is 

.2 



while 



d 



•\ 2 



^ - \l + -K{z)\\ogi 



(56) 
(57) 
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Going over to cosmic time ^ = t°(^)+2, gives the metric in Kasner form with exponents 

P.(.) = ^^. (59) 
The following relations are immediately verified: 

/D-l 



D-l \ 
Y,Pa{z) -1] logt, 
a=l J 



(60) 



D-l 



E pI (^) = 1- (61) 



The behaviour of the fields near the singularity is thus of Kasner type, modified, as usual, 
by the presence of the dilaton. Note that, at the two tips of the singularity, and k diverge in 
such a way that the Kasner exponents near the tips are simply pi = 1, pi = 0. This corresponds 
to a (contracting) Milne-like metric which, being non-singular, nicely matches the non-singular 
behaviour in regions II and III. Away from these two points, k and can take any value: it 
is easy to verify that the whole Kasner sphere can be covered by appropriately choosing the 
initial data. 

For generic z, the collision leads to a singularity showing all the characteristics of a cos- 
mological singularity. The Kasner exponents depend on z and therefore the behaviour of the 
metric will depend on the point of the singular hypersurface that we approach as ^ — > 0. It 
is therefore interesting to study the signs of the Kasner exponents for different choices of the 
coefficients and n. In particular, we would like to see whether, far from the tips of the 
singularity, inflation may take place. 

By the definition of a{z), we see that the denominators of the Kasner exponents ( |5BD and 
(|59D are always positive. Concerning the numerators, we note that a {z) is a quadratic form in 
ej and k. The equation a{z) =0 thus defines an ellipsoid in the space (e^; k) with pi positive 
outside of this ellipsoid and negative inside. For each i the equation Pi = defines the plane 
K + ei = passing through the centre of the former ellipsoid. The p^'s are positive on the 
semispace containing the positive k semiaxis. Finally, we have to remember the constraint on 
the sum of the ipi, which becomes X^^i = in the asymptotic regime. 

In conclusion, we see that there indeed exists a region where all the Kasner exponents are 
negative. In fact, if we stay inside the ellipsoid on the side where k < 0, and take all close 
enough to zero, all the exponents are negative and we have infiation in all directions. In the 
opposite situation (i.e. outside the ellipsoid, with k > and close to zero) all exponents are 
positive and we have contraction in all directions. All intermediate cases are also possible as 
we can let any number of be large and negative and take at least one large enough and 
positive to balance the others. 
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6 Colliding waves with antisymmetric field 



Although we have written the equations in a manifestly 0((i, (i)-covariant form, we have only 
managed, so far, to find solutions in region IV for vanishing antisymmetric field. We can 
perform, however, 0{d,d) boosts on our solutions to introduce it. This procedure can be 
worked out in any number of the dimensions. Our purpose in this section is to discuss the 
main features of the cosmologies obtained by this procedure and to clarify the role of the 
antisymmetric field, particularly in connection with the work of Ref . |^ . 



Let us write the transformation matrix Q in terms of four d- dimensional blocks 



P R 
Q S 



The constraint (11) reads 



P^Q + Q^P = 0, 
R^Q + = Id, 
P^'S + Q^'R = Id, 
R^S + S^R = 0. 



From eq. (0), with M in the form (44), we obtain the inverse of the new metric G' 

G'-^ = P^G^^P + Q^GQ 
and extract also the antisymmetric field 

B' = (r^G-^P + S^GS) G'. 



(62) 



(63) 
(64) 
(65) 
(66) 



(67) 



(6^ 



It is interesting to investigate the asymptotic behaviour of the new cosmologies obtained 
this way. We try to be very general, considering an initial metric that is diagonal and has the 
Kasner form 

G = diag(t^i,...,t^<^), (69) 

with the Aj ordered so that i < j implies that |Aj| > |Aj|. In section 5 we showed that the 
asjTuptotic approach to the singularity of the colliding waves metric is of this type, if we let 
the Aj depend on z. 

When we perform the O {d, d) transformation, the new metric G' is no longer diagonal. 
However, we can diagonalize it by a step-by-step procedure. Exploiting the fact that G has the 
form (^), the eigenvalue equation for G'^^ takes the following asaymptotic form 

E E {P^ar^'Pib + Q^at^'Q^h) Vb = aVa- (70) 

b i 
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Assuming, in a first instance, that all Aj are non-vanishing, we define two new matrices U 
and U', whose components are 




Pij ifA, >0 
Qij if Ai < 



Qij if Xi > 
Pii if Ai < 



(71) 



(72) 



The dominant term in the eigenvalue equation as t ^ becomes 



(73) 



b 



This is the eigenvalue equation for the projector to Ui, i.e. the first eigenvalue of G'^^ is 
a = |f/ipt~l'*'^l with eigenvector Va = Uia- Since the vectors orthogonal to Ui annihilate the 
dominant term we just considered, in order to find the other eigenvalues and the corresponding 
eigenvectors, we consider the successive sub-dominant terms in the original eigenvalue equation 
in the space orthogonal to Ui. 

The next term is 



The eigenvector of this term is v = U2- however, in order to annihilate the dominant term, we 
have to take the component of U2 (asymptotically) orthogonal to Ui as the second eigenvector. 

The corresponding eigenvalue is a = (^|f/2p — -^^j^j^it-) t~^^^^. The next eigenvector must now be 
taken in the space orthogonal to f/i and U2, and so on. In this way, we can build the complete 
system of eigenvectors of G'~^, combining the rows of U, which are rows of the two matrices P 
and Q. The d eigenvalues are then proportional to t"'^*' and therefore they are all diverging as 
t — >■ 0. It is easy to extend this discussion to the case of some vanishing A,. This clearly implies 
asymptotically vanishing eigenvalues for G. 

In order to complete our discussion, we have to mention some possible exceptions to the 
above-mentioned situation. Not always, the rows of U, used to build the eigenvectors, form an 
independent set of vectors. For a particular class of 0{d, d) transformations, having det f/ = 0, 
the result is different. It is easy to understand the modifications: each time we face a Ui linearly 
dependent on the previous rows, we just ignore it and proceed with the algorithm to the next 
term in G'~^. When we have exhausted the first d dominant terms, we will still have to find 
one eigenvalue (having skipped one term in the diagonalization) . Taking the next term of G'~^, 



we can find the missing eigenvalue. In conclusion, we can see that, if the rank of [/ is s < c?, 
then G'"^^ will have s diverging and d — s vanishing eigenvalues. Note that, even if the matrix P 
is chosen freely, the matrix Q is related to P by eq. (|63|). We have to check whether, because of 
this relation, the determinant of U is necessarily vanishing for some initial values of the Kasner 
exponents. 



^U2at ^^^^U2bVb = aVa. 



(74) 



b 



in the form 




(75) 
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First, we observe that (^) implies 

detPdetg = (-l)'^detQdetP. (76) 

Therefore, for odd d, either the determinant of P or the determinant of Q vanish. We may 
now distinguish two classes of 0{d, d) transformations. The identity transformation belongs to 
those giving a vanishing det Q, and the same is true for all transformations that can be reduced 
to the identity in a continuous way. For metrics that are inflating in all directions before 
the transformation, U = Q and the transformations having det Q = 0, instead of turning all 
eigenvalues to contracting ones, leave at least an inflating. Conversely, for metrics contracting 
in all dimensions, U = P and transformations with det P = cannot leave all the eigenvalues 
unchanged, but induce inflation in at least one direction. 

Consider, as the last peculiar case, the one in which the exponents in the initial metric are all 
positive except one. Then it is possible to show that the constraint forces det U = if P is 
non-singular. The same argument can be repeated for metrics inflating in all directions except 
one. In this case, U is composed by — 1 rows from Q and one from P. If Q is non-singular, 
then det [/ = and again we have one inflating dimension left. 

Summing up, we can say that cosmologies with antisymmetric field that can be obtained, 
by an 0{d,d) transformation, from a metric having Kasner behaviour near the singularity, 
generally contract in all dimensions. In even dimensions, at least one inflates when the original 
metric has just one inflating or just one contracting dimension. In odd dimensions, if the 0{d, d) 
transformation has det Q = 0, then, starting from a full inflating metric or from a metric with 
one inflating dimension, we are left with one inflating dimension. If det P = 0, then one 
dimension inflates when we start from a full contracting metric or from a metric having just 
one contracting dimension. 

We should also recall that none of the 0{d, d) transformations affects the gu component of 
the metric, which, therefore, can either inflate or contract. 



The results of Ref. |]T6| are compatible with ours, since, for d = 2, our statements, when the 



determinant of G inflates, can be summarized as a change of sign in A while ^p is left invariant. 



7 Conclusions 

In this paper, extending previous work we have modelled the onset of pre-big bang inflation, 
from asymptotically trivial initial conditions as the result of the collision of two plane waves, 
made of gravitons, dilatons, and Kalb-Ramond massless particles, in any number D = d + 2 of 
space-time dimensions. We showed that the evolution of the system is described in terms of a 
compact and elegant set of 0{d, (i)-invariant equations and that properties of the solutions can 
be studied in full generality in three of the four regions defined by the planar-collision problem. 
This already enables us to argue that the formation of a curvature singularity in the future (to 
be identified with the big bang) is generic. 
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However, so far, we have not been able to solve the general problem analytically in the 
fourth, and most interesting, region, except for the case of dilatonic and parallel-polarized 
gravitational waves. Nonetheless, using 0{d, d) transformations, we were able to construct new 
solutions containing the KR field and to discuss the physically relevant properties of these new 
solutions. The main conclusions that appear to emerge are the following: 



While the formation of a singularity is generic, the existence of inflating regions near the 
singular surface is only generic in the absence of the KR form. 

The KR form, at least when generated from 0{d,d) transformations, tends to generate 
contraction rather than expansion, in agreement with other results and arguments . 

Since our equations appear to reduce, near the singularity, to those studied previously 



in the homogeneous (or quasi-homogeneous) case, it looks very likely that the 



above behaviour will persist in the general solution. 



Eventually, one may be led to the conclusion that the most generic APT initial conditions 
(i.e. those containing all possible kinds of waves in the initial state) can hardly produce such 
a flat, homogeneous and isotropic Universe to dispense us completely from the more standard 
kind of potential-energy-driven post-big bang inflation. Actually, it is very likely that, after 
exit from pre-big bang inflation, the dilaton and other moduli will find themselves dispaced 
from the minima of their non-perturbatively-generated potentials and that further inflation 
will result from their rolling down towards them. A similar conclusion seems to follow from 
completely different, more phenomenological arguments, i.e. from the recent data analysis of 
CMB anisotropies at small angular scales |T^, which appears to confirm the need for adiabatic 
perturbations of the kind naturally provided by potential-driven inflation, but absent in the 
PBB scenario. 

If so, we will have to back up from the early claims that the PBB scenario can replace 
altogether standard inflation, and settle instead for a complementary role it would play in 
providing the initial conditions that standard inflation badly needs, and in "explaining" , from 
more natural and generic initial conditions, the most mysterious event in the entire life of our 
Universe, the big bang. 

Note Added 

While this paper was being written, we received a new paper by Damour and Henneaux 



1^] containing a discussion of the outcome of 0{d, d) transformations on Kasner-like solutions. 
Their results, obtained with a different diagonalization procedure, agree with those discussed 
in our section 6. 



Acknowledgements 

We are grateful to Maurizio Gasperini for help during the early stage of this work and to 
Thibault Damour for discussions concerning his recent work, Refs. P, |18|. V.B. was supported 
by fund ex 60% D.P.R. 382/80, MURST, Italy. He also wishes to thank the Theory Division 
at CERN for hospitality while this work was carried out. 



14 



References 



[1] J.E. Lidsey, D. Wands and E.J. Copeland, [hep-th/9909061| , to appear in Phys. Rep.; G. 



Veneziano, |hep-th/0002094| , and references therein 



[2] M. Turner and E. Weinberg, Phys. Rev. D56 (1997) 4604; 

N. Kaloper, A. Linde and R. Bousso, Phys. Rev. D59 (1999) 043508; 

A. Buonanno, K.A. Meissner, C. UngareUi and G. Veneziano, Phys. Rev. D57 (1998) 2543; 
M. Gasperini, |hep-th/0004Ti9| . 



[3] A. Buonanno, T. Damour and G. Veneziano, NucL Phys. B543 (1999) 275. 

[4] D. Christodoulou, Comm. Math. Phys. 105 (1986) 337, Comm. Pure AppL Math. 54 (1991) 
339 and Comm. Pure AppL Math. 56 (1993) 1131. 



[5] A. Feinstein, K.E. Kunze, M.A. Vazquez-Mozo, |hep-th/0002070| . 
[6] T. Damour and M. Henneaux, |hep-th/0003139| . 

[7] V.A. Behnskii, I.M. Khalatnikov and E.M. Lifshitz, Adv. Phys. 19 (1970) 525. 
[8] P. Szekeres, J. Math. Phys. 13 (1972) 286. 

[9] K.A. Meissner and G. Veneziano, Phys. Lett. B267 (1991) 33 and Mod. Phys. Lett. A6 
(1991) 37. 

[10] A. Sen, Phys. Lett. B271 (1991) 295; 

S.F. Hassan and A. Sen, NucL Phys. B375 (1992) 103. 

[11] M. Gasperini and G. Veneziano, Phys. Lett. B277 (1992) 256. 

[12] G. Veneziano, Phys. Lett. B265 (1991) 287; 
A.A. Tseythn, Mod. Phys. Lett. A6 (1991) 1721; 
A. A. Tseythn and C. Vafa, NucL Phys. B372 (1992) 443. 

[13] J. Maharana, J.H. Schwarz, NucL Phys. B390 (1992) 3. 

[14] A. Buonanno et aL, ReL 0. 

[15] U. Yurtsever, Phys. Rev. D37 (1988) 2803. 

[16] D. Clancy, A. Feinstein, J.E. Lidsey, R. Tavakol, Phys. Rev. D60 (1999) 043503. 



[17] A.E. Lange et aL, Mro-ph/000"5004 



A. Balbi et aL, |astro-ph/OOlM^. 



[18] T. Damour and M. Henneaux, [hep-th/ 0006 171 . 



15 



